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From the measurements of F. E. Giesecke on the flow of 
liquids through submerged orifices two definite laws have 
been deduced. When the diaphragm is completely closed, 
the liquid flows backwards along the wall of the tube, and 
forward in a central vein, the radius of which has been 
calculated for laminar and turbulent flow. The theory by 


v. Karman for turbulent flow has been applied to a tube of 
circular cross section. It is shown that both the backward 
flow and the jet issuing from the orifice are unstable. By an 
experiment it is made evident that a jet may sometimes 
describe a spiral motion. 





ROFESSOR F. E. Giesecke! has made very 

careful measurements of the flow of water 
through submerged orifices, from which definite 
laws can be deduced. He measured the loss of 
head as a function of the velocity and of the 
diaphragm, shown in Fig. 1. The diameter of the 
pipe was 2.65 cm, the diameter of the orifice 
varied from 0.75 to 2 cm, and the velocity of the 
water varied between 2 to 200 cm per second. 
The results were plotted in logarithmic coordin- 
ates. Using, however, ordinary coordinates, 
drawing as abscissae the square of the velocities, 
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1F. E. Giesecke, ‘‘Loss of Head in Submerged Orifices,”’ 
Heating, Piping and Air Conditioning, Vol. 2 (1930), 
p. 611. 


v®, and as ordinates the loss of head Ah, I find the 
straight lines of Fig. 2, according to which 


Ah= Cv’, (1) 


a law which follows immediately from the prin- 
ciple of dimensional homogeneity. The straight 
lines of Fig. 2 correspond to the various diameters 
of the diaphragm, from 0.4 to 0.8 inch. 

If I plot on logarithmic coordinate paper the 
loss of head Af as function of the diameter, then 
I find the curves of Fig. 3, which are nearly 
straight lines showing, however, a slight curva- 
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ture for small values of Ah. These curves repre- 
sent Ah for various velocities. The straight line 
part of these curves can be represented by 
log Ah= —ax+b or (x = D=diameter); or 


Ah = Ahye~*”. (2) 
Combining Eqs. (1) and (2) we obtain 
Ah=Cv%e-, (3) 


where C is a constant and a=tan a is the slope of 
the curves of Fig. 3. The average value of a is 
73° 15’, hence a=tan a= 3.32. When the straight 
lines of Fig. 3 are prolonged until they intersect 
the axis of ordinates, we find the points corre- 
sponding to the loss of pressure when the dia- 
phragm is closed; it is easily found that Ah/v!-76 
is constant. But 1.75=7:4 and _ Prandtl’s 
theory of turbulent motion gives for the shearing 
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force zo at the wall of the tube: 
to= Cpvtr-ty7/4, 


If the diaphragm is completely closed, then the 
pressure increases and a pressure wave travels 
through the liquid column, with the velocity of 
sound; but there is also another phenomenon 
recently observed by E. E. Ambrosius and J. C. 
Reed, namely, the liquid returns at the dia- 
phragm and moves along the wall of the tube 
in a direction opposite to the original motion in 
the center. Let in Fig. 4 the radius of the tube be 
R, the radius of the central bore be 7, its average 
velocity be 7, and the average velocity over the 
tubular cross section 7(R?—7r,?) be a, then when 
the diaphragm is closed, we must have in every 
cross section 


tri; =7(R?—1r;") ie, (4) 


expressing the fact that the resultant flow 
through any cross section is zero. 

If we assume at first laminar flow, we have fs 
the steady state of flow 


pate —(R?—1). (5) 
l An 


The maximum velocity u, in the center of the 
circular cross section is equal to 


_pi —p2 1 |. 
l An 





Um 


while the average velocity 7% is equal to 
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gion ee (6) 
lL 8m 2 





The profile of the motion is a parabola, shown in 
Fig. 5. For a certain radius 1 the velocity is 
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equal to the average velocity. This radius is de- 
termined by (5) and (6) as 


R?—r?=R?/2, or 17;2=R?/2. (7) 
Using (6) and (7) we obtain the difference 


Unt’? — Un (R?—r,?) = 2aar?—ain(r;?) =aar,’. 


“The volume of the cylinder (Fig. 4) (ABCD) 


minus the volume of the ring (AA’, BB’, EE’, 
FF’) is equal to the volume of the cylinder ABEF, 
but this is also equal to ABCD—EFCD; hence 
the volume of EFCD is equal to the volume of 
AA’, BB’, EE’, FF". lf we therefore assume that 
i over the ring AA’ BB’ is equal and opposite to 
a over the circular cross section 7r,2, then the 
average resultant velocity through the whole 
cross section is zero; in other words, the average 
velocity %, must be equal to #.= 4 in formula (4) 
which then reduces to an identity. We have 
therefore immediately after closing the dia- 
phragm 


r1=R/v2=R-0.71. (8) 


‘Hence the mass which flows through 7(R?—7r,?) 


backwards, is equal to the mass which flows for- 
ward through zr’; hence also the average mo- 
mentum and the average energy through the two 
cross sections per unit time is the same. 

But the average velocity 7 cannot remain con- 
stant, because the energy will gradually and 
even rapidly become dissipated in the form of 
heat. Indeed, the work done by the forces of fric- 
tion per unit time and unit volume will be equal 
to 


dW /dt = idp/dx = n8ri?/x R?, 


and per unit length equal to 7R°adp/dx = n8rit” 
and the work done through dx per unit time will 
be dW = 8rni*dx, or dW /dx = 82rni? is equal to the 
decrease of the kinetic energy per unit time of the 
volume 7R?-1, which gives 


82rnii? = — (d/dt)(prR?u?) ; 
hence 
a= f2e~16nt/ PR? (9) 
We can also find the decrease of the average 
velocity along the axis or the tube in the following 


way. The difference in kinetic energy entering 
and leaving the volume dx7R? per unit time is 





equal to the kinetic energy in this volume trans- 
formed in heat. The entering kinetic energy is 
equal to 7R?(p/2)#?a; the kinetic energy leaving 
is equal to tR*(p/2)(a#+0a0x/dx)?(i+dtdx/dx); 
the difference is equal to (— 37R®pu?/2)(du/dx)dx, 
and this is transformed into heat. [8xna#dx], 
hence #=—16nx/3pR?+i%. We shall see 
moreover that this sort of laminar motion is 
unstable. But before considering that problem, 
we shall consider the steady turbulent flow in a 
tube of circular cross section with radius R. If 
we apply the theory of turbulent flow of Th. v. 
Karman? to a tube, then we obtain the formula 


1 
t= tint -(ro/0)!Clog (1 (r/R))+(7/ "J, (10) 


where k is a constant between 0.36 and 0.4; p is 
the density, 7) the shearing stress at the surface; 
the average velocity then becomes equal to 


1 R 2 
-—|{ u2ardr=— { urdr 
aR? 0 R? 


2 21 
-— f unrdr += (10/9) f [log (1—(r/R)}) 
+(r/R)*|rdr. 


Integrating by expansion in a series, I obtain 
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2 Th. v. Karman, Nachrichten d. Ges. d. Wiss. Géttingen 
(Jan., 1930), p. 58. 
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This average value @ is equal to u at a distance "; 
from the center, where 7, is determined by the 
equation 


wnemnsMToe(-G))*C) 


77 1 To 4 
242) 
60k\ p 


or by the equation 


77 


Zt (1-(3) (5) 1 


Solving this transcendental equation for 7, I find 
by trial and error 


r1=R-0.766, (12) 


instead of the value for ~ for laminar motion 
given in Eq. (8). The profiles of the laminar and 
the turbulent motion are given in Figs. 5 and 6 
‘for maximum velocities of 64 and 720, respec- 
tively. Extrapolated from the tables of I. 
Nikuradse,* ro p= 965. 

If a flowing column of water is suddenly 
stopped, it has been assumed that the whole 
kinetic energy per unit volume 1/2 m7? is trans- 
formed into elastic energy of compression, the 
pressure increase Ap at the diaphragm being 
given by 


Ap=AV-K=p—0=p), 


where K is the coefficient of compressibility, and 
AV the change of the volume. The conservation 
of energy gives: 
tmi?=pAV, but p=p/2 
=>): p/K=p?/2K 
or 
3]. Nikuradse, Forschungsarbeiten 289, 38 (1927). 
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(13) 


K } 
p=a(—m-m) =tum(V,)}, 
m 


where V,=(K/m)! is the velocity of sound in 
the column, # is the sudden pressure increase at 
the diaphragm. Now the question arises, whether 
this pressure increase is the same, when the liquid 
turns around and moves backwards. We have 
seen that the kinetic energy, flowing through any 
cross section 7R? is equal to zero, so that indeed 
the whole original kinetic energy before closing 
the valve is transformed again into elastic energy 
of compression according to Eq. (13) though from 
an analogy with a solid elastic body, rebounding 
from a rigid wall it might appear at a first thought 
that the pressure p might be twice as large. But 
the analogy is not complete, because in the case 
of the moving liquid column only one-half of the 
water rebounds. 

We have considered one liquid column flowing 
inside another one in opposite direction. Such a 
motion is unstable. This follows from Rayleigh’s* 
criterion of stability of jets, applied to our pres- 
ent problem. We consider two coaxial cylinders 
of liquids, (Fig. 7) the outer one having the 
density p’, the velocity V’ and the velocity poten- 
tial 


¢’=C’ cosh k(r—re)ei(ott*), 
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while for the inner cylinder 
¢=C cosh k(r+rj)eottk, 


and the disturbance at the boundary h is equal to 
h=He'**)), k=2x/d, w=2nv, where v is the 
frequency, A the wave-length; then the motion is 
stable if 


(g(p—p') +k*T)(p cot hkri+p’ cot hkre) 
—kpp’ cot hkr, cot hkr.( V—V")?>0. 


T is the surface tension between the two non- 
viscous liquids, g the acceleration of gravity. In 
the present case p=p’, T=0, V=V’, hence the 
condition of stability cannot be satisfied for any 
value of \. Even if V were not equal to — V’, and 
if \ were very small, still no stability would ob- 
tain. This conclusion also holds for the jet, issuing 
from the submerged orifice. These phenomena 
have recently been studied by Ambrosius and 
Reed. Behind the orifice there appears for small 
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velocities a system of small vortices accompany- 
ing the central jet, shown in Fig. 1b. For higher 
velocities the central jet begins to oscillate, 
showing a system of vortices, Fig. 1a, which re- 
minds me somewhat of the vortex ally of v. 
Karman. But it may be that there are vortex 
rings surrounding the central jet, vortex rings 
whose axes are inclined toward the direction x 
of the general motion. These vortices disappear 
gradually as the liquid flows downward, giving 
rise to steady turbulent motion. It is also possible 
to imagine that the central jet describes a spiral 
motion. In order to test this possibility, an ex- 
periment was made which is illustrated by Fig. 8. 
A test tube is filled with turpentine; a wick, 
passing through a capillary tube, burns with a 
small steady flame in the lower end of a glass tube 
of 5 cm diameter and 30 cm length. A very nar- 
row column of smoke rises from the flame, often 
showing a spiral motion. Sometimes the jet forms 
a ribbon, which appears twisted like an elastic 


band. 
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Viscous Properties of Polyisobutylene 


Joun D. Ferry* anp GeorGE S. Parks, Department of Chemistry, Stanford University 
-(Received July 22, 1935) 


The viscosity of polyisobutylene (molecular weight about 
4900) has been studied from 70°C to 175°C with a falling 
sphere viscometer and from —53°C to 24°C with a con- 
centric cylinder viscometer. In the measurements with the 
latter apparatus, elastico-viscous effects appeared, and the 
viscous behavior was non-Newtonian. The apparent vis- 
cosities were accordingly extrapolated to zero torque to 
give values for the viscosity at zero shear. The viscosity 


ranges from 19 poises at 175°C to 10” poises at —53°C, 
showing a rather low temperature coefficient. At — 76°C, 
the midpoint of the “transition region” for polyisobutylene 
glass, the viscosity has an extrapolated value of about 10% 
poises, in agreement with a general rule for glass-forming 
materials. Macroscopic and microscopic interpretations of 
viscosity anomalies are discussed. 





INTRODUCTION 


LY MERIC organic compounds whose mole- 

cules are supposed to have the form of 
prolonged chains of repeated units have been 
the subject of much experimental study in dilute 
solution. Among the properties of these linear 
macromolecules which have been investigated are 
diffusion,! streaming double refraction,? sedi- 
mentation,’ and dipole moments.‘ Above all, the 
effect of long, rod-shaped molecules on the vis- 
cosity of a solvent has been extensively studied 
experimentally®: © and has been given consider- 
able theoretical treatment.’ On the other hand, 
the physical properties of the pure polymers in 
the absence of solvent have received little 
attention. 

In connection with an investigation® of the 
physical properties of a sample of polyiso- 
butylene, with special attention to the formation 
of a glass at low temperatures, we have studied 
the viscosity of this substance in the pure form 
over a range of temperature from —53°C to 
175°C. Measurements were made by the falling 
sphere method and by the Couette or concentric 


* Shell Research Fellow at Stanford University for the 
academic year 1934-35. 

1R. O. Herzog and H. Kudar, Zeits. f. physik. Chemie 
A167, 343 (1934). 

2 R. Signer and H. Gross, Zeits. f. physik. Chemie A165, 
161 (1933). 

3E. O. Kraemer and W. D. Lansing, J. Am. Chem. Soc. 
55, 4319 (1933); J. Phys. Chem. 39, 153 (1935). 

4W. Gallay, Koll. Zeits. 57, 1 (1931). 

5H. Staudinger, Die Hochpolymeren Organischen Ver- 
bindungen (Julius Springer, Berlin, 1932). 

6 E. O. Kraemer and F. J. Van Natta, J. Phys. Chem. 
36, 3175 (1932). 

7E, Guth and H. Mark, Ergeb. d. exakt. Naturwiss. 12, 
115 (1933). 

8]. D. Ferry and G. S. Parks, to be published in the 
J. Chem. Phys. 


cylinder method. These viscosity determinations 
were undertaken for comparison with previous 
measurements on other glass-forming liquids®: 
of simple molecular structure and low molecular 
weight. 

The present study also provides examples of 
anomalous or ‘“‘non-Newtonian”’’ viscous be- 
havior. In the temperature range in which the 
concentric cylinder method was applied, the 
apparent viscosity for a given temperature and 
torque was a function of the time (due to elastic 
effects), approaching a constant value (cf. Filon 
and Jessop," Stott,!? and Umstatter'*) ; whereas 
the final, constant value was, at a given tempera- 
ture, a function of the torque applied (cf. Reiner," 
Lawrence,'® and Philippoff'*). 


MATERIAL 


The sample of polyisobutylene used in this 
work had an average molecular weight of 4900 
as determined by cryoscopic measurements in 
benzene solution; actually the sample was a 
mixture of polymers with a wide range of 
molecular weights. 

The molecules of polyisobutylene, like those of 
other linear polymers, are generally supposed to 
be unbranched chains of great length. This 
hypothesis, while not at present susceptible to 


°G. S. Parks, L. Barton, M. E. Spaght and J. W. 
Richardson, Physics 5, 193 (1934). 

10(G. S. Parks and M. E. Spaght, Physics 6, 69 (1935). 

uL. N. G. Filon and H. T. Jessop, Phil. Trans. Roy. 
Soc. A223, 89 (1922). 

12. H. Stott in W. S. Turner, The Constituticn of Glass 
(Society of Glass Technology, Sheffield, 1927), pp. 73, 155. 

13H. Umstitter, Koll. Zeits. 70, 174 (1935). 

14M. Reiner, Physics 5, 321 (1934). 

15 A, C. S. Lawrence, Proc. Roy. Soc. A148, 59 (1935). 

16 W. v. Philippoff, Koll. Zeits. 71, 1 (1935). 
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proof, is supported by the above-mentioned 
studies in dilute solution,': *: 5 ° especially those 
of Staudinger. The most probable arrangement 
for the isobutylene units in polyisobutylene is 


~ CH; CH; 
| | 
ial Ga iil 
| | 
CH; CH; CH; 


If bond lengths and angles are the same as in 
normal paraffins, a molecular weight of 4900 
corresponds to a length of about 220A. However, 
more complicated arrangements may also be 
involved, with side chains longer than the methyl 
groups. 

The density of our sample, determined at 
three temperatures, is given in Table I. 


TABLE I, Density of polyisobutylene. 


Temperature 23.4°C 100.2 178.2 
Density 0.9074 g/cm’ 0.873 0.835 


MEASUREMENTS BY THE FALLING SPHERE 
METHOD 


Viscosity determinations over the temperature 
range from 70°C to 175°C were carried out by 
the falling sphere method, with the same 
apparatus previously employed for measuring 
the viscosity of glucose. The formula of Laden- 
burg,” in which the effect of the sides and ends 
of the cylindrical column of liquid is taken into 
account, was used to evaluate the viscosity from 
the measured rates of fall of steel bearing-balls. 

Under the circumstances of these experiments, 
there was no evidence for non-Newtonian viscous 
behavior. Comparative determinations at the 
same temperature with different sizes of balls 
agreed within experimental error ; and successive 
determinations at the same temperature showed 
no trend indicative of irreversible changes or 
persistent anisotropic arrangements induced by 
the motion of the spheres. It was concluded, 
therefore, that a Newtonian flow obtained, and 
that the quantities calculated by application of 
the Ladenburg formula represented ‘‘true’’ vis- 
cosities. 

Over seventy individual determinations were 
made in the given temperature range. In a plot 


17 R, Ladenburg, Ann. d. Physik (4) 23, 447 (1907). 
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of the logarithm of viscosity against temperature, 
most of the points fell within a distance from the 
smooth curve which corresponded to a deviation 
in the viscosity of 4 percent. Measurements made 
before and after heating to the highest tempera- 
ture employed (175°C) were in satisfactory 
agreement, showing that no decomposition or 
oxidation sufficient to affect the viscosity was 
occasioned by the heating. Duplicate determina- 
tions were also made at 21°C and were in 
reasonable agreement with the results from the 
concentric cylinder measurements described be- 
low. Representative values of the viscosity and 
its logarithm are given in Table II. These values 


TABLE II. Viscosity of polyisobutylene by the falling sphere 
method. 








VISCOSITY, 7 VISCOSITY, 7 








Temp. (°C) —POISES LoGw 7 | Temp. (°C) POrIsES LoGw 7 
21 5.8 104 4.76 120 1.26 x 10? 2.10 
70 1,63 x 10% 3.21 130 85 1.93 
80 9.1 «10? 2.96 140 59 1.77 
90 5.1 X10? 2.71 150 42 1.62 
100 3.1 10? 2.49 160 30 1.48 
110 1.9 x10? 2.28 170 22 1.35 








are represented in Fig. 4 by the unbroken portion 
of the curve, the scale being too small here to 
show the individual determinations by the falling 
sphere method ; except those at 21°C, which are 
represented by a triangular point. 


MEASUREMENTS BY THE CONCENTRIC CYLINDER 
METHOD 


No determinations of viscosity were made 
between 24°C and 70°C. Below 24°C, measure- 
ments were carried out by the Couette or con- 
centric cylinder method. 


Apparatus 


The apparatus employed was very similar to 
that previously used in the studies on glucose,® 
but modified to permit measurements at low 
temperatures. A Pyrex test-tube (internal di- 
ameter 2.27 cm) and a hollow steel rod (diameter 
0.95 cm) constituted the outer and inner cyl- 
inders, respectively ; the latter was immersed in 
the polyisobutylene to a depth of about 7 cm. 
Temperature control was accomplished by a 
double jacket consisting of two cylindrical copper 
buckets separated by an air gap of 1 cm. The 
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inner jacket, which was equipped with an 
adequate stirrer and a heating coil, contained a 
bath of ethyl alcohol surrounding the Pyrex 
test-tube. The outer jacket was suspended in a 
3-liter Dewar jar, and was packed with ice or 
with carbon dioxide snow. The whole arrange- 
ment was rigidly attached to a small rod of 
channel iron, the bearings of the rotor or inner 
cylinder being accurately centered. To minimize 
leakage of heat down the channel iron into the 
bath, the apparatus was enclosed in a cabinet of 
Celotex. In order to attain the lowest tempera- 
tures, small quantities of liquid air were intro- 
duced into the space between the two copper 
jackets. The lowest temperature at which a 
viscosity measurement was made was —53°C. 
The lowest temperature at which satisfactory 
thermal equilibrium could be maintained for a 
period of an hour or two (long enough for a com- 
plete series of determinations for a consistency 
curve) was —43°C. The temperature of the 
bottom of the steel rotor was measured by a 
copper-constantan thermocouple in connection 
with a White potentiometer. A similar thermo- 
couple measured the temperature of the alcohol 
bath in the inner jacket. 

Torque was applied to the inner cylinder by 
exerting a force at the rim of a wheel of known 
dimensions. This force was supplied by a weight 
attached to a thread which was wound around 
the wheel and passed over a pulley. The problem 
of the frictional drag of the rotor bearings was 
considered carefully. The apparatus was cali- 
brated when empty by weighting equally both 
ends of a thread wrapped around the wheel, and 
determining the additional weight required to 
be added to one end to overcome frictional 
resistance to turning. The frictional moment was 
found to be separable into two terms, the first 
independent of load (attributed to the vertical 
reaction on the bearings) and the second pro- 
portional to the total load and amounting to 
6 percent of it (attributed to the horizontal 
reaction on the bearings). In all viscosity de- 
terminations, the data were corrected by sub- 
tracting from the actual moment applied the 
corresponding frictional moment. 

Other details of construction and operation of 
the viscometer were similar to those in the 
glucose study previously mentioned. The usual 
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formula employed for calculation of the vis- 
cosity is 
n=kC/w, (1) 


where C is the applied torque, w the angular 
velocity of the rotor, and k is a constant cal- 
culable from the dimensions of the apparatus.® 
This equation depends on the assumption that 
Newtonian flow obtains. In the case of non- 
Newtonian flow, as in the present experiments, 
application of the formula gives “apparent 
viscosity”’ values. The latter probably do not 
describe the fundamental viscous properties of 
the material, but show only how it behaves in 
the given type of experiment. A more funda- 
mental quantity is the viscosity at zero shear, 
obtained by extrapolation of the apparent vis- 
cosity, as explained below. 


Dependence of deformation rate upon time: 
Elastico-viscous effects 


Under a given load, the rate of turning was a 
maximum when torque was first applied to the 
rotor; it fell off with time to approach a steady 
value. This behavior was undoubtedly due to 
combined elastic. and viscous displacements 
caused by the shearing stresses in the polyiso- 
butylene (elastico-viscous effect'®). When the 
torque was removed, the rotor ‘‘unwound”’ to a 


‘certain extent; the rate of this elastic return, 


also, was rapid at first and fell off with time 
asymptotically to zero. 
In a few cases, the course of the initial elastic 
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Fic. 1. Angular displacement of inner cylinder plotted 
against time, showing elastic effects. Portion A of the 
curve represents elastic and viscous displacement ; portion 
B, viscous displacement alone; and portion C, elastic 
return. 














deformation when the torque was applied and 
that of the final elastic return when the torque 
was removed were followed in detail.!* Curves of 
angular displacement against time, plotted from 
these data, were all of the type shown in Fig. 1. 
Each curve consisted of a portion A representing 
combined elastic and viscous deformation, a 
portion B representing viscous deformation 
alone, and a portion C representing elastic 
return. In the majority of cases, however, no 
data were taken for the portions A and C of 
such curves. The load was simply applied for a 
time long enough to satisfy the elastic displace- 
ment, i.e., until successive determinations of the 
angular velocity of rotation were in agreement. 
At the lowest temperatures, it was necessary to 
wait as much as twenty minutes before the rate 
of rotation attained this stationary value (repre- 
sented by the slope of portion B of the displace- 
ment-time curve). 


Dependence of elastic effects upon torque applied 


Extrapolation of the straight portion B of the 
displacement-time curve to zero time gives an 
angle @;, which represents the excess initial 
displacement over what would have occurred had 
the viscous flow of portion B been the only type 
of deformation from the start. This angle, 0:, and 
that of the total elastic return upon removal of 
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Fic. 2. Elastic displacement (initial ‘and final) in 


radians, plotted against torque in dyne-cm. Temperature 
= 242.5°K. 


1® We are indebted to Mr. B. L. Crawford, Jr., for 
assistance in some of these measurements. 
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Fic. 3. Consistency curve for polyisobutylene: Angular 
velocity in radians per second plotted against applied 
torque in dyne-cm. Temperature = 257.8°K. 


load, 42, are related to the elasticity of the 
material. The fact that the two are not found to 
be identical suggests that the elasticity of the 
flowing material is different from that of the ma- 
terial at rest. The dependence of @; and 62 upon 
the applied torque is illustrated in Fig. 2, which 
gives data from a single series of experiments. 
The initial displacement, 6;, is roughly pro- 
portional to the torque (corresponding to Hooke’s 
law) ; the final elastic return, 62, is much smaller 
than 6, at low torques. Both quantities appear to 
fall off from the proportionality law at high 
torques, suggesting that some sort of elastic 
limit has been reached. 


Dependence of apparent viscosity upon torque 
applied: Non-Newtonian flow 


At a given temperature, the stationary value 
of angular velocity of turning was not pro- 
portional to the torque applied. Curves of angular 
velocity against torque were all of the type 
shown in Fig. 3, which represents a ‘‘consistency 
curve,” in the terminology of Reiner,!* and is a 
good example of non-Newtonian flow. In this 
diagram, apparent viscosities (nc’) are repre- 
sented by the reciprocal slopes of lines drawn 
from the origin to the consistency curve (multi- 
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19 Reiner’s ‘‘consistency curve’ (reference 14) is actually 
a plot of rate of shear against shearing stress; in the 
Couette apparatus, these variables are proportional, 
respectively, to the angular velocity of turning and to the 
applied torque. 
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Fic. 4. The logarithm of the viscosity of polyisobutylene 
as a function of temperature. The circular points represent 
measurements by the concentric cylinder method; the 
unbroken line above 340°K and the triangular point, 
measurements by the falling sphere method. 


plied by the apparatus constant, k). The vis- 
cosity of the material at zero shear is given by 
the reciprocal slope of the tangent to the con- 
sistency curve at the origin, multiplied by R. 

The non-Newtonian flow of polyisobutylene is 
in marked contrast to the strictly Newtonian 
flow observed for other glass-forming liquids. 
whose consistency curves are straight lines 
passing through the origin.®: '° 

The viscosity at zero shear is defined as 


no=lim kRC/w=lim ne’ (2) 
C0 c+0 


and is a property of the material at rest, in the 
absence of abnormal molecular arrangements 
induced by viscosimetric experiments. Hence 
this is the quantity to be compared with the vis- 
cosities measured at higher temperatures by the 
falling sphere method (where the shear was in- 
sufficient to occasion non-Newtonian effects), 
and with viscosities of Newtonian systems in 
general. Accordingly, at a number of tempera- 
tures, data for consistency curves were obtained, 
and values for no were calculated by extrapolating 
nc’ to zero torque. The ratio nc’ /no varied from 
0.9 to 0.8 for the smallest torque used, and was 
as low as 0.3 for the greatest torque. The ratio 
nc’/no for a given torque varied less than 20 
percent over a temperature range from 230°K to 
270°K, in which no varied two thousand-fold. 
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At temperatures where no consistency curves 
were determined, the apparent viscosity nc’ for 
some given torque was corrected to give no by a 
coefficient estimated from preceding and suc- 
ceeding determinations. 


Representative values for viscosity at zero shear 


The final values of the logarithm of 0 for 
twenty-three viscosity determinations by the 
concentric cylinder method, with two different 
portions of the polyisobutylene sample, are 
plotted in Fig. 4. Most of the points fall within a 
distance of the smooth curve corresponding to a 
deviation in the viscosity of less than 12 percent. 
Representative values of the zero-shear viscosity 
and its logarithm are given in Table III. 


TABLE III. Values of no for polyisobutylene by the concentric 
cylinder method. 











Temp. (°C) no LoGuw 0 | Temp. (°C) no LoGio no 
— 50 6.5 < 10° 9.81 —10 2.5108 6.40 
—40 6.2 x 108 8.79 0 6.9 10° 5.84 
— 30 7.6107 7.88 10 2.3105 5.36 
— 20 1.3107 7.10 20 8x 104 4.90 








THE VISCOSITY-TEMPERATURE CURVE 


The viscosity data obtained from both meth- 
ods: of measurement are plotted against the 
absolute temperature in Fig. 4, where the un- 
broken line above 340°K and the triangular 
point show the results of the falling sphere 
method, and the circles represent the values of 
no from the concentric cylinder experiments. 
This curve has the same form as observed for 
other glasses,® 1° except that it is much less 
steep. The temperature coefficient of the viscosity 
of polyisobutylene is rather small; the quantity 
d\n mo/d(1/T) is for polyisobutylene on the 
average only 0.3 that for glucose, and 0.4 that 
for boron trioxide (comparisons made at corre- 
sponding viscosities). 

In the previous studies on glucose, boron 
trioxide, and glycerol®® (in the last two cases by 
extrapolation from the measurements to lower 
temperatures), the viscosity of each of these 
glass-forming liquids was found to have a value 


of between 10% and 10 poises at the mid- 


20G. S. Parks and W. A. Gilkey, J. Phys. Chem. 33, 
1428 (1929). 
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point of the “transition region’’ of temperature 
where the heat capacity and thermal expansion 
coefficient undergo sudden increases with rising 
temperature. This rule appears to hold irre- 
spective of the nature of the material or of the 
mean temperature at which the transition takes 
place (513°K for boron trioxide; 187°K for 
glycerol). Measurements of the heat capacity 
and thermal expansion of polyisobutylene® place 
the midpoint of the transition interval for this 
material at 197°K (i.e., —76°C). Extrapolation 
of the viscosity curve to this temperature gives 
a value of 10%? poises. This conformity of 
polyisobutylene to the general rule is note- 
worthy, in view of its great molecular size and 
anisodimensionality (in contrast to the other 
glass-forming liquids), and its exhibition of non- 
Newtonian viscous behavior. 


DISCUSSION 


The anomalous viscous phenomena exhibited 
by polyisobutylene have been observed for many 
other systems, and numerous formulations and 
interpretations have been put forward to de- 
scribe them, without significant success. These 
theoretical treatments fall into two classes, the 
macroscopic, which consider only macroscopically 
measurable properties of the material in question, 
and the microscopic, which involve speculations 
concerning the microscopic structure of the 
material. 

The starting-point of many of the macroscopic 
theories is the fundamental equation of Maxwell, 


ds/dt=(1/E)dF/dt+F/n, (3) 


where s is the deformation, F the force pro- 
ducing it, 7 the viscosity and E the elasticity. 
Under conditions where the two terms of the 
right member are of the same order of magnitude, 
“elastico-viscous” effects of the type observed 
(Fig. 1) are to be expected. They have been 
found for pitch,” glucose glass fibers,” silicate 
glasses,’ and celluloid." Quantitative inappli- 
cability of the Maxwell equation is, however, 
shown by the fact that it predicts an exponential 
decay of elastic strains by viscous flow, whereas 





*1F, T. Trouton and E. S. Andrews, Proc. Phys. Soc. 
London 19, 47 (1909). 


22 Unpublished experiments by Mr. J. D. Reagh in this 
laboratory. 
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the observed decay does not follow an ex- 
ponential law.?* Umstitter," starting from Eq. 
(3), deduced a completely opposite elastic effect 
from that observed for the above-mentioned 
substances; it should result in a rate of deforma- 
tion which is low at first and rises to approach a 
steady value. Michelson* noted that the decrease 
in the rate of elastic deformation with time did 
not follow an ordinary exponential law, but that 
it could be fitted to an exponential law involving 
the square root of the time. The data for polyiso- 
butylene fit this formulation fairly well, but its 
physical significance is obscure. Filon and 
Jessop" employed a formulation involving 2’, 
but it is purely empirical, and it does not give a 
deformation rate which approaches a constant 
value. Thompson*® has derived an elaborate 
theory of visco-elasticity, satisfactorily tested by 
experiments on rubber. Unfortunately none of 
his formulations is applicable to our type of ex- 
periment (deformation under constant torque). 
Eisenschitz?® deduced from (3) that whenever 
(F/E)? was not negligible compared with unity, 
non-Newtonian flow should result, and the 
following expression should hold approximately : 
nr’ /no=1/(1+ F?/E?*), (4) 
where 7,’ is the apparent viscosity for a force F. 
This formulation does not hold for the data on 
polyisobutylene, whose consistency curves (Fig. 
3) cannot be fitted to an equation of the form 
w=aC+bC (the equivalent of Eq. (4)). Besides, 
the ratio n/n has been found to vary but little 
with temperature, demanding by (4) a similar 
constancy in E; whereas as a matter of fact the 
latter varies considerably with temperature. 
Philippoff'* modified the equation of Eisenschitz 
and introduced a new parameter, 7,,, the viscosity 
at infinite shear. However, the data for polyiso- 
butylene cannot be well fitted to his formulation 
either. 
The microscopic explanations for the origin of 
non-Newtonian flow have been summarized by 
Reiner.'* Of his long list, the most important in 


23E. Hatschek, The Viscosity aw Liquids (G. Bell and 
Sons, London, 1928), pp. 224-227. 

4A A, Michelson, J. Geol. 25, ‘405 (1917). 

% J. H.C. Thompson, Phil. Trans. Roy. Soc. A231, 339 
(1933); J. C. Eccles and J. H. C. Thompson, Proc. Roy. 
Soc. A148, 171 (1935). 

26 R, Eisenschitz, Koll. Zeits. 64, 184 (1933). 
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the case of polyisobutylene are probably orienta- 
tion and stretching out of tangled macro- 
molecules. Alignment of molecules in the lines 
of flow, either by orientation or by uncoiling and 
untangling, would diminish the viscous resistance 
to flow; and both these effects would increase 
with the shear, giving rise to a non-Newtonian 
consistency curve. More specifically, the orienta- 
tion would depend on the rate of shear, and the 
stretching on the shearing stress. On the other 
hand, it seems probable that the elastic behavior 
is due to the stretching phenomena alone, and 
hence that the elastic properties and the viscous 
anomalies of the material do not have an 
identical microscopic origin, and cannot be con- 
nected so simply as Eisenschitz postulates. 

The explanation of elastico-viscous effects for 
polyisobutylene cannot depend solely on its 
molecular size and anisodimensionality, because 
these effects occur in materials of simple mo- 


lecular structure, such as glucose. However, a 


much greater shearing force is required to evoke 
elastico-viscous effects in glucose than in poly- 
isobutylene; the elastic coefficient of polyiso- 
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butylene is rather smaller than those of other 
glass-forming liquids, compared at corresponding 
viscosities. It is possible that, at high shearing 
stresses, glucose would exhibit non-Newtonian 
flow; that is, the E/isenschitz criterion for 
existence of non-Newtonian behavior may be 
qualitatively applicable. 

A complete explanation of viscosity anomalies 
must consider Brownian movement, which com- 
petes with the orienting and stretching forces 
and renders them less effective at higher tem- 
peratures. The absence of non-Newtonian effects 
in the viscosity experiments on polyisobutylene 
by the falling sphere method is attributed to the 
considerable thermal agitation at those tempera- 
tures and to the low shearing stresses involved 
in the flow. 

It would be desirable, in anticipation of an 
adequate theory of viscous behavior, to make a 
complete experimental study of the dependence 
of elastic and viscous flow on time, shearing 
stress, and temperature, using preferably a one- 
component system, of which many glass-forming 
liquids offer convenient examples. 
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Compressional Waves in Media with Complex Viscosity 


A. GEMANT, Engineering Laboratory, Oxford University 
(Received July 16, 1935) 


The author recently suggested the introduction of a 
complex viscosity in the theory of vibrations of materials 
which lie between the purely crystalline and the liquid 
state. This conception is utilized here in the theory of com- 
pressional waves in an infinite medium. One of the conclu- 
sions, namely, the decrease of the attenuation factor with 
increasing frequency, is in agreement with observations. 


Another conclusion is an increase of the propagation veloc- 
ity with frequency. Both changes occur in the neighborhood 
of a certain characteristic frequency, given by the reciprocal 
of the so-called relaxation constant. Below this frequency 
range the material vibrates like a liquid; above this range 
like a crystal. 





HE behavior of a medium with regard to 

longitudinal compressional waves is gen- 
erally described by means of two constants. 
One of them is the velocity of propagation, c, 
the other the attenuation factor, a. Apart from 
a term in the latter caused by heat conduction, 
a is proportional to the viscosity of the medium, 
u. This viscosity can be calculated from the 
attenuation factor, and the agreement between 
this value of the viscosity and its value known 
from other determinations is usually satisfactory, 
at least in order of magnitude. Such is the case 
with measurements concerning gases and liquids. 
Usually these media are dealt with since they 
allow the conditions corresponding to propaga- 
tion of waves in an infinite medium to be 
reproduced satisfactorily. Even more frequently 
they are investigated in tubes, since the theory 
of waves in tubes is well developed. 

More difficult is the investigation of solids, as 
their finite shape, that is the consideration of the 
boundary conditions, makes the evaluation of 
measurements rather complicated; correspond- 
ingly, only a small amount of data is available. 
It seems, however, that such measurements show 
a certain discrepancy between theory and experi- 
ment. In the first place, the derived value of the 
viscosity is much lower in order of magnitude 
than the static value. Secondly, the viscosity 
seems to decrease with increasing frequency of 
the waves. 

For example, J. H. C. Thompson! compares 
the values for u obtained by several authors for 
metals and glass, as follows: Hettwer? obtained 
with very slow torsional vibrations the order 


1 J. H. C. Thompson, Phil. Trans. Roy. Soc. A231, 339 
(1933). 


2 Hettwer, Wien Ber. 134, 51 (1925). 


10'5, an order probably not far from the static 
value. Honda and Konno,’ investigating the 
torsional pendulum at higher frequencies, found 
the order 10%. S. L. Quimby‘ applied ultrasonic 
vibrations set by a quartz oscillator to solid rods 
and derived for uw the order 10*. The decrease of 
the order of the viscosity with increasing fre- 
quency is very definite and especially the last 
order obtained by Quimby is much lower than 
the static value. 

Crandall’ mentions in his Vibrating Systems 
and Sound unpublished experiments of P. W. 
Bridgman and H. M. Trueblood, who investi- 
gated the propagation velocity and damping of 
sound waves in rubber. They also observed 
qualitatively a behavior similar to that described 
above. 

In a recent paper, J. C. Eccles and J. H. C. 
Thompson® measured the damping of longi- 
tudinal vibrations in rubber for frequencies of 
the order 10%. The value for yu as calculated for 
pure black rubber ranges from 8.10* for a 
frequency 16, to 5.10° for a frequency 32. The 
order for yw is again certainly lower than that of 
the static value. The order 3.10* is that of the 
static viscosity of a rather viscous oil, as can be 
seen from measurements of the author.’ This 
value was observed for some lubricating oils at 
0°C; the static value for rubber must be con- 
siderably higher. 

These peculiarities can be explained by a 

8’ Honda and Konno, Phil. Mag. 42, 115 (1921). 

4S. L. Quimby, Phys. Rev. 25, 559 (1925). 

5]. B. Crandall, Vibrating Systems and Sound (Van 
Nostrand Co., 1927), p. 111. 

6 J. C. Eccles and J. H. C. Thompson, Proc. Roy. Soc. 
A148, 171 (1935). 


7A. Gemant, 


Viskositat u. Fliessfestigkeit Berl. Ber. 
(1932). 
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theory suggested recently by the author.* ° 
The idea of this theory was given by considera- 
tions of a somewhat analogous peculiarity of the 
viscosity in the field of certain organic colloids, 
generally called elasto-viscous liquids. In order 
to explain their behavior, the Maxwellian con- 
ception of elasto-viscosity has been applied to 
them by H. Hencky!® and others with consider- 
able success. The theory of the author is a 
generalization of this idea. The generalization is 
twofold. Instead of assuming a viscous resistance 
and a resilience connected in parallel, as in the 
original conception, it is assumed that the 
damping term is the resultant of a more compli- 
cated combination of units, and that it consists 
in the simplest case of three elements at least. 
Such units are known from the treatment of 
plastic bodies.!! This combination will be called 
the complex viscosity and denoted by V. 

The second generalization refers to the appli- 
cation. It can be applied to any kind of vibration 
in viscous media. Thus the application to the 
flow of liquids, as has been done by Hencky and 
others, is only one possibility. Another applica- 
tion is concerned with electric losses in amor- 
phous solid dielectrics, giving an extension of 
the theory of Debye” to such materials. This 
extension has been carried out quantitatively in 
one of the above-quoted papers.° 

A third application is concerned with the 
acoustical waves in solids, the attenuation of 
which cannot be explained on the basis of a 
simple viscosity constant yu, as has just been 
pointed out. 

In this note only the simplest case, really the 
fundamental one, is discussed, namely, the 
damping of waves in a laterally unlimited iso- 
tropic medium. Another important case, namely, 
waves in a medium contained in a rigid circular 
tube, will be treated in another communication. 
Again, the case of vibrations in a solid bar needs 
separate treatment. 

In order to fix our ideas, it should be recalled 
that V is considered in the simplest case as a 
parallel combination of a viscous resistance 7, 

8 A. Gemant, Naturwiss. 23, 406 (1935). 

*A. Gemant, Conception of a Complex Viscosity and its 
A pplication to Dielectrics, Trans. Faraday Soc., in press. 

10H. Hencky, Ann. d. Physik 2, 617 (1929). 

First Report on Viscosity and Plasticity, Acad. of 


Sciences at Amsterdam, 1935. 
2% P, Debye, Polar Molecules (New York, 1929). 
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on the one hand, and of a series combination of 
a resilience E and another viscous resistance 7, 
on the other. Here 7 is the static viscosity and 
corresponds physically to an agitational dis- 
placement of molecules, such as is generally 
assumed to take place in liquids. E is the 
elastic resistance to a regular displacement of 
molecules without disturbance of their structure, 
a motion which characterizes displacements in 
crystals. 7 is the frictional resistance to be over- 
come during this crystalline displacement. In a 
so-called Newtonian liquid 7 is so small that the 
other branch in V becomes unimportant; in a 
crystal, on the other hand, 7 is so large that 
only the branch E—7 is of importance. The 
general complex impedance covers therefore the 
whole transition state between a pure liquid and 
a pure crystal. To this state belong all the 
amorphous or glassy materials, and all those 
which can be considered (in a certain temperature 
range) as a colloidal mixture of crystalline and 
liquid or glassy constituents (paraffin wax, for 
instance). 


The quantity To= 10o/E (1) 


is called the visco-elastic time constant of 
relaxation for reasons shown below. It is easy 
to show® that V may be expressed as a complex 
quantity: 


V=7'—-iE’/w, (2) 


where w is 27 times the frequency. Both 7’ and 
E’ are dependent on frequency. Considering 
that in all practical cases 


no>; (3) 
we get for n’ with good approximation: 


n’ = nolL 1+(n/E) row ]/[1+(row)?] = (4) 
and for E’: E’=E(row)?/[1+(7w)?]. (5) 


It may be seen that for frequencies below 
1/70, n’ becomes m and E’ very small, whereas 
for frequencies above 1/70, 7’ becomes and 
E’ : E. This is the reason why 7» is called the 
time constant of relaxation; the transition from 
an original elastic strain to a viscous flow takes 
place in an interval of the order 7». 

We have to apply Eq. (2) to the theory of 
propagation of longitudinal waves. Denoting the 




















displacement in the direction x by &, the velocity 
by &, we have: 


ae 40°F AE 
———e—— c— =0, (6) 
0 3 0x dx? 


where p=density and x=bulk modulus of 
elasticity. 

(6) is the classical equation derived by H. 
Lamb. Since the quantity V, which has to 
replace the simple viscosity u is a complex one, 
Eq. (6) has to be brought to a complex form. 


This is done by assuming for the solution: 
= Eetott—slen, (7) 


where c;=complex propagation constant. 
In substituting (7) in (6), we obtain 


—p+i(4/3)pw/cPr+x/cr=0. (8) 
Now u has to be replaced by V, and we get: 
— p+1(4/3)n'w/c?+(4/3)E'/cr+x/cr=0, 
cP =(x/p)(1+(4/3)E’/x+1(4/3)n'w/x). 


If the damping term is small compared with 
the others, we have to a first approximation : 


kK 4 E’\ 7} 
ONO 
p 3 kK 
The complex expression for c; permits calculation 


of the real propagation velocity and the attenu- 
ation factor. Since (7) is identical with 


hence 


13.9’ w/k 


). (9) 
1+(4/3)E’/« 





£= Eye azeivlt—2/0) (7a) 
we have the identity: 
€,=c(1+iac/w). (10) 
Comparing (10) with (9), we have for c, 
c=[(x/p)(1+(4/3)E’/«) }}, (11) 
and for a: a= 3n'w/ pc, (12) 


where E’ and 7’ are given by (4) and (5). 
The following conclusion can be drawn from 
these equations. For w<1/79: E’<E and 7’ =». 
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Since E is of similar order to x, this means that 
c=(x/p)' and a=%3nw?/pc*. These equations are 
the same as obtained by considering the medium 
as a liquid with bulk elasticity x and viscosity mo. 
Below the critical frequency the medium behaves 
like a liquid. 

For frequencies w>1/70, we have E’=E and 
n’ =n. Hence 


c=[(x/p)(1+(4/3)E/x) }', 


and a= 3nw*/ pc. 

The medium here behaves like a crystal. The 
bulk elasticity is now (x+4£/3), a quantity 
which ‘includes the elasticity owing to lattice 
distortion, as defined by (2). The damping is 
given by the dissipation accompanying crystal 
distortion. If E is small compared with x, the 
change in c may be only a second-order effect. 

The transition between the two cases just 
considered occurs continually in the frequency 
range~ 1/79. Now in a number of cases, especially 
with solids, such as glass at room temperature, 
To is certainly larger than 1 sec. In other cases, 
like rubber, it might be less than 1 sec. For 
liquids showing structure viscosity it is certainly 
considerably smaller than 1 sec. The range in 
which the greatest variation of the propagation 
velocity (i.e., the most pronounced dispersion of 
waves) and of the attenuation factor is to be 
expected depends on the value of 7». 

Whereas the theory seems to explain the 
above-mentioned peculiarities satisfactorily, a 
quantitative check would require more experi- 
mental data. It should be added only that the 
experimental curve relating the variation of 
damping with frequency is probably more flat 
than indicated by (4). This would mean that 
the simple equivalent circuit as accepted for V 
is not adequate. There are experimental indica- 
tions for the exponent of the argument (w7 9) in 
Eqs. (4) and (5) being nearer to unity. 

From the standpoint of comparison with 
experiment, the extension of the theory to 
propagation in tubes and to vibration of bars is 
of great importance. 





